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Current quantum computing architectures lack the size and fidelity required for universal fault- 
tolerant operation, limiting the practical implementation of key quantum algorithms to all but the 
smallest problem sizes. In this work we propose an alternative method for general-purpose quantum 
computation that is ideally suited for such “prethreshold” superconducting hardware. Computa¬ 
tions are performed in the n-dimensional single-excitation subspace (SES) of a system of n tunably 
coupled superconducting qubits. The approach is not scalable, but allows many operations in the 
unitary group SU(n) to be implemented by a single application of the Hamiltonian, bypassing the 
need to decompose a desired unitary into elementary gates. This feature makes large, nontrivial 
quantum computations possible within the available coherence time. We show how to use a pro¬ 
grammable SES chip to perform fast amplitude amplification and phase estimation, two versatile 
quantum subalgorithms. We also show that an SES processor is well suited for Hamiltonian simu¬ 
lation, specifically simulation of the Schrodinger equation with a real but otherwise arbitrary nxn 
Hamiltonian matrix. We discuss the utility and practicality of such a universal quantum simulator, 
and propose its application to the study of realistic atomic and molecular collisions. 

PACS numbers: 03.67.Lx, 85.25.Cp 


I. INTRODUCTION AND MOTIVATION 
A. The promise of quantum computation 

A universal quantum computer, if one could be built, 
would transform information technology by providing 
vastly increased computational power for certain special¬ 
ized tasks, such as quantum simulation [IHi and prime 
factorization Superconducting electrical circuits 

operating in the quantum regime [^, 0 have emerged as 
an extremely promising platform for realizing a large- 
scale, practical machine. Yet the quantum algorithms 
actually demonstrated to date—with any architecture— 
have been limited to only tiny, few-qubit instances El¬ 
ls^. In superconducting circuit or circuit QED imple¬ 
mentations, which benefit from the inherent scalability 
of modern solid-state electronics, this barrier in qubit 
number does not reflect any limitation of the underly¬ 
ing device fabrication or infrastructure requirements, but 
rather that larger problem sizes would also require longer 
computations and hence additional coherence. Quantum 
algorithms typically have (uncompiled) circuits that are 
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spatially narrow but temporally very deep. In this work 
we propose an alternative approach to superconducting 
quantum information processing that allows one to cir¬ 
cumvent this restriction and realize much larger compu¬ 
tations within the available coherence time. 

A general-purpose quantum computer that is useful for 
practical applications must, of course, be error corrected 
and scalable. The standard model of an error-corrected 
quantum computer is the gate-based fault-tolerant uni¬ 
versal quantum computer, where “errors” acting on all 
device components and at any step during the computa¬ 
tion can be corrected as long as they are weak enough— 
below an error threshold J3j-l3g —and not highly corre¬ 
lated in space or time [37h4(]|| . Scalability means that the 
number of physical qubits required to perform a partic¬ 
ular computation—the physical volume of the quantum 
computer—scales as a polynomial function (preferably 
linear) of the problem size. It also means that it is actu¬ 
ally possible, in practice, to add more qubits. 

A realistic picture of an error-corrected superconduct¬ 
ing quantum computer based on the surface code [4l|, 
is beginning to emerge [i^ . The surface code is the most 
practical, best performing fault-tolerant approach known 
to date, and is especially amenable to implementaion 
with superconducting circuit technology. However, the 
resources required for a practical machine are consider- 
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able: Fowler et al. [i^ estimated that factoring a 2000-bit 
number would require about 2x10® physical qubits, using 
Beauregard’s modular exponentiation 0 and a surface 
code quantum computer operating at 99.9% fidelity. If 
there was no decoherence or noise, and no errors of any 
kind to correct, then it would be possible to factor an 
A^-bit number with the Beauregard algorithm using only 
2N + 3 ideal qubits, or 4003 ideal qubits in the case con¬ 
sidered. Thus, error correction imposes a physical qubit 
overhead of 2 x 10®/4003 Ri 5 x 10"'^. Note that in quan¬ 
tifying the error-correction overhead here we distinguish 
between ideal (error free) qubits—the fictional entities 
usually appearing in quantum algorithms—and logical 
qubits, which must also include the many additional an- 
cillas necessary for fault-tolerant gate implementation. 
Similarly, You et al. [i^ estimated that it would take 
about 5 X 10® physical qubits to calculate the ground 
state energy of a 100-spin transverse-field Ising model to 
99% accuracy using the same 99.9%-fidelity surface code 
quantum computer. This well known statistical mechan¬ 
ics model maps especially well to a quantum computer, 
and for N spins on a line would require only -|- 1 ideal 
qubits for a calculation of the ground state energy (using 
iterative phase estimation). So the physical/ideal ratio in 
this quantum simulation example is 5xl0®/101 « 5x10'^, 
the same as for factoring. Therefore we expect that, 
in practice, surface code error correction will impose an 
overhead of 

# physical qubits _ . 

# ideal qubits 

where we have allowed for future optimization and other 
improvements. Crudely, a factor of about 10 in the over¬ 
head estimate comes from replacing ideal qubits with 
enough logical qubits to both encode those ideal qubits 
and to distill the auxiliary states needed to perform fault- 
tolerant operations on them, and a factor of about 10® 
comes from replacing each logical qubit with enough 
physical qubits to enable a sufficiently long computation. 

B. Prethreshold quantum computation 

The complexity of building even a small fault-tolerant 
universal quantum computer suggests that this objective 
may take some time to achieve. In the meantime, ex¬ 
perimental quantum information processing is limited to 
either the very small problem sizes discussed above, or 
to nonuniversal approaches such as analog quantum sim¬ 
ulation 11^-0 , quantum annealing [0 , or other 

special-purpose methods [5l|, [0. In this work we la¬ 
bel any quantum computation without an error-corrected 
universal quantum computer as prethreshold, referring to 
the threhold theorems of fault-tolerant quantum compu¬ 
tation, because exceeding a fidelity threshold is a neces¬ 
sary condition for large-scale error correction. 

Table U compares three broad approaches to quantum 
computation with prethreshold hardware: Small system 


refers to gate-based computations with a few qubits, 
which have been used to test fundamental concepts of 
quantum information processing, demonstrate hardware 
functionality, and assess qubit and gate performance. 
The SES method is also general purpose, but should en¬ 
able quantum speedup (this is discussed below). However 
neither approach is scalable. Analog quantum simulation 
and other scalable, special-purpose approaches trade uni¬ 
versality for a faster route to speedup. 

The SES quantum computer described in this work is 
universal in the sense that it can implement any gate- 
based algorithm or quantum circuit. As a simulator it 
can directly emulate any (real) Hamiltonian, including 
time-dependent Hamiltonians. When we refer to a simu¬ 
lated Hamiltonian in this context we mean a Hamiltonian 
written in some basis—a real, symmetric matrix with no 
special structure. We assume that the Hamiltonian ma¬ 
trix to be simulated has been specified externally, as is 
typically the case when using a classical computer. The 
SES processor solves the Schrodinger equation defined by 
this Hamiltonian. 


TABLE I. Three approaches to prethreshold quantum com¬ 
putation and simulation. The left column lists the attributes 
achievable by an error-corrected universal quantum computer. 



small system 

SES method 

analog QS/spec purp 

scalable 

X 

X 

% 

universal 

% 

% 

X 

speedup 

X 

% 

% 

arb accuracy 

X 

X 

X 

arb runtime 

X 

X 

X 


The superconducting SES method introduced here has 
features in common with the siMle-photon protocols of 
Reck et al. 0 and Cerf et al. [0, as both use only one 
excitation, and therefore do not utilize genuine entangle¬ 
ment. The optical realization uses a recursive algorithm 
to first decompose a given nx n unitary U of interest 
into a sequence of SU(2) beam-splitter transformations. 
This decomposition determines an arrangement of beam 
splitters, phase shifters, and mirrors, that will unitar- 
ily transform n input ports (optical modes) to n output 
ports according to the desired U. However, the supercon¬ 
ducting realization is better suited for quantum simula¬ 
tion than the optical approach because the Hamiltonian 
is directly programmed. In particular, to optically sim¬ 
ulate Schrodinger evolution under a given Hamiltonian 
matrix H, one would have to first compute the evolution 
operator on a classical computer, and then decom¬ 

pose it into beam-splitter transformations, but avoiding 
the classical computation of (or the time-ordered 

exponential if H is time dependent) is the motivation for 
the quantum computation in the first place [0. Nei¬ 
ther the superconducting SES nor the single-photon op¬ 
tical approaches are scalable—they both require expo¬ 
nential physical resources—and should not be considered 
as viable alternatives to the standard paradigm of error- 
corrected universal quantum computation. But they are 
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both suitable prethreshold methods. 

II. QUANTUM COMPUTATION IN THE SES 

A. Hardware model: The programmable SES chip 

Consider the following model of an array of n coupled 
superconducting qubits, 

= + ® ’ ( 2 ) 
i ii' 

written in the basis of uncoupled-qubit eigenstates. Here 
i,i' = 1,2,..., n, and 



The Ci are qubit transition energies and the gu' are qubit- 
qubit interaction strengths; both are assumed to be tun¬ 
able. (Factors of h are suppressed throughout this pa¬ 
per.) gai is a real, symmetric matrix with vanishing di¬ 
agonal elements. We also require microwave pulse control 
of at least one qubit, and simultaneous readout (projec¬ 
tive measurement in the diagonal basis) of every qubit. 
The model m describes a fully connected network or 
complete graph of qubits, which we refer to as an SES 
processor. This should be contrasted with local quan¬ 
tum computer models that have coupling only between 
nearby qubits (nearest neighbors, for example). The SES 
method can be applied with a wide variety of qubit- 
qubit interaction types (see Appendix but without 
loss of generality we restrict ourselves here to the simple 
coupling of ([U. Alternatively, tunably coupled 


B. Single-excitation subspace 

The idea we explore in this paper is to perform a quan¬ 
tum computation in the n-dimensional single-excitation 
subspace of the full 2"-dimensional Hilbert space. This is 
the subspace spanned by the computational basis states 

|i) =cJ|00---0) = |0---l,---0), (4) 

with i = 1, 2,..., n. We call the set of |i) the SES basis 
states. It is simple to prepare the quantum computer in 
an SES basis state from the ground state |00 • • • 0), and 
it will remain there with high probability if the following 
conditions are satisfied: 

1. The coupling strengths \gu'\ are much smaller than 
the Ci, which is usually well satisfied in supercon¬ 
ducting circuits. 

2. Single-qubit operations such as tt and 7r/2 rotations 
about the x or y axes are not used during the com¬ 


resonators (with tunable frequencies) can be used instead 
of qubits [5^ . Although we assume an architecture based 
on superconducting circuits (or circuit QED), the SES 
method might apply to other future architectures as well. 

The quantum computer model © might be consid¬ 
ered unscalable, because of the 0{n^) tunable coupling 
circuits and wires, a position that we also adopt here. 
In gate-based universal quantum computation, the fully 
connected and local quantum computer models are equiv¬ 
alent in the sense that any quantum circuit implemented 
by a fully connected quantum computer can be imple¬ 
mented by a local quantum computer after adding chains 
of SWAP gates, which only introduce polynomial over¬ 
head. However, this equivalence is restricted to the stan¬ 
dard gate-based approach and does not apply here. 

Superconducting qubits have been reviewed in Refs. Q 
and Q. Although the model ([2]) can be realized with 
several qubit designs, the transmon qubit currently 
has the best performance For concreteness we 

assume a qubit frequency e/27r in the range of 5.45 to 
5.55 GHz and coupling strength g/2Tr in the range —50 
to 50 MHz. An n-qubit SES processor also requires 
n(n — l)/2 coupler circuits and the associated wires or 
waveguides. A variety of tunable couplers can be used 
for this purpose. Here we consider a modification of the 
tunable inductive coupler developed by Chen et al. [6l| 
for superconducting Xmon qubits; this design has been 
demonstrated to implement tunability without compro¬ 
mising high coherence. Our modification replaces the 
direct electrical connection of each qubit to a coupler cir¬ 
cuit wire with an inductive coupling to the wire, which 
scales better. An SES chip layout that avoids excessive 
crossovers is illustrated in Fig.[I] The tunable interaction 
strength g for this coupler design is derived in Appendix 

m 


putation. However, 27r rotations are permitted and 
are very useful (these can be implemented as z ro¬ 
tations, which do not require microwaves), tt rota¬ 
tions about X or y can be used to prepare SES basis 
states from the system ground state |00 • • • 0). 


3. The quantum computation time is significantly 
shorter than the single-qubit population relaxation 
time Ti. 


An SES pure state is of the form 


'0) = ^a*|z), ^|aip = l. 

2=1 2=1 


( 5 ) 
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FIG. 1. (Color online) Circuit layout for n = 5 SES processor, with the crosses representing Josephson junctions. Each 
horizontal circuit is an Xmon qubit with capacitance C, tunable junction inductance Lj, and n — 1 additional coils (each with 
self-inductance Lq and mutual inductance m) for coupling to other qubits. Dotted lines indicate dc and microwave control 
lines for each qubit, as well as readout circuits. Each coupler wire contains a Josephson junction with inductance Lc tuned by 
a magnetic flux $. Control lines for SES matrix elements gi 2 , ■ ■ ■ ,gi 5 are also indicated. This circuit is discussed further in 
Appendix iBl 


For example, the states ([S]) include the maximally entan¬ 
gled VF-type state 


unif) 


|1) + |2) + ■ ■ ■ + |n) 

y/n 

| 10 ... 0 ) + | 01 ... 0 ) + ...+ 100 ... 1 ) 

y/n 


Although the state ([6]) is entangled, and could be used 
to violate Bell’s inequality, the entanglement is somewhat 
artificial as there is only one “particle”. 


C. SES Hamiltonian 

The advantage of working in the SES can be under¬ 
stood from the following expression for the SES matrix 
elements of model ([2]), namely 

'^ii' — “t“ Qii'• (7) 

Because the diagonal and off-diagonal elements are di¬ 
rectly and independently controlled by the qubit frequen¬ 
cies and coupling strengths, respectively, we have a high 
degree of programmability of the SES component of the 
quantum computer’s Hamiltonian. This property allows 
many n-dimensional unitary operations to be carried out 
in a single step, bypassing the need to decompose into ele- 
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mentary gates, and also enables the direct quantum sim¬ 
ulation of real but otherwise arbitrary time-dependent 
Hamiltonians. However, we have some restrictions: 

1. T-Lii' is always real, whereas the most general Hamil¬ 
tonian matrix is complex Hermitian. The exper¬ 
imentally available control parameters, consisting 
of n qubit frequencies and n{n — l)/2 coupling 
strengths, are sufficient to control the n{n -|- l)/2 
independent parameters of an n x n real symmetric 
matrix. 

2. There are experimental limitations on the range of 
values that the and gar can take. We define t^max 
to be the magnitude of the largest coupling avail¬ 
able in a particular experimental realization; a cur¬ 
rent realistic value is about 50 MHz. 

We will leave the discussion of possible generalizations to 
complex SES Hamiltonians for future work. The limita¬ 
tions on the ranges of the and gui do not, by them¬ 
selves, limit the class of Hamiltonians that can be simu¬ 
lated, because a model Hamiltonian intended for simula¬ 
tion is first rescaled to conform to that of the SES chip 
(this is explained below). 

It will be useful to refer to a “typical” SES Hamiltonian 
"H, which we assume to have the following properties: 
H is a real, symmetric matrix with each element taking 
values in the range —5max to ^max- This form follows from 
© after removing an unimportant term proportional to 
the identity matrix, 

n^n- ( 8 ) 

where Wref is a convenient (possibly time-dependent) ref¬ 
erence frequency. Then we assume that the qubit fre¬ 
quencies Ci can be tuned within i^max of tOrei, and we 
assume that the couplers can be tuned between —^max 
to (7max- A possible choice for Wref is the mean qubit 
frequency (l/n)X)iei- 

It will also be useful to consider the statistical prop¬ 
erties of an ensemble of typical SES matrices. We can 
always write a time-independent SES Hamiltonian in the 
standard form 

T-L = 5max K, (9) 

where K is real symmetric matrix with every element 
satisfying 

-1 < Kw < 1 . ( 10 ) 

We define a real random matrix ensemble of dimension n 
as follows: The n diagonal elements are independent ran¬ 
dom variables Ku, each uniformily distributed between 
— 1 and 1. The n(n—1)/2 elements Ki^i/ are independent 
random variables also uniformily distributed between — 1 
and 1. The remaining elements Ki^ii are fixed by the 
symmetry requirement. The standard deviation of each 
element is <7^ = l/-\/3. 



matrix dimension n 


FIG. 2. (color online) Spectral bandwidth of K matrices ver¬ 
sus n. Data (open circles) are averaged over 1000 random 
instances of K. The solid line is the function 1.58xn° ®®. The 
dashed line is the function (A/y/S) x^/n. 



matrix dimension n 


FIG. 3. (color online) Spectral bandwidth for larger n. The 
solid line is the function 2.06 xn°'®^. The dashed line is the 
function {A/ y/i)Xy/n. 

The first property we study is the mean spectral band¬ 
width of itr, the difference between the largest and small¬ 
est eigenvalues. Let Ai < A 2 < ■ • • < A„ be the ordered 
eigenvalues of K. From the Wigner semicircle law 
we expect that, in the large n limit, 

An - Ai-)> = ^n5, (11) 

where the overbar denotes averaging over the ensemble 
defined above. The bandwidth of typical SES states in 
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an n-qubit processor therefore scales at large n as 

-^max -^min ^ 5max n 2 . (12) 

In Figs. [H and [3] we plot the simulated bandwidth of K 
as a function of n, and compare the simulation data with 
the asymptotic form (HU). From Fig. [2] we conclude that 
for modest SES matrix sizes, 

i^max - E^in « 1-58 g^ax (13) 

The second property we study is the mean level spacing 
of K. Let 

n—1 

AA =-- Ai+i — Xi (14) 

n — 1 

be the mean spacing between adjacent eigenvalues. Av¬ 
eraging d over the ensemble defined above, we expect 
that in the large n limit 

(15) 

n V J 

In Fig. 0] we plot the simulated average level spacing of 
K as a function of n, and compare the simulation data 
with the asymptotic form m- From Fig. |4]we conclude 
that for modest SES matrix sizes, 

A^«1.89g„,axn-°-4®. (16) 

The results (fTSl) and (fTHl) give two relevant energy scales 
present in a typical SES spectrum. 

Any unitary quantum circuit or operation acting on q 
qubits can be mapped to and implemented on an SES 
chip with n = 2^ qubits (this exponential growth of n 
is what makes the SES method unscalable). We can 
say that the SES processor simulates the g-qubit sys¬ 
tem, with the advantage of being able to perform multi¬ 
qubit operations in a single step. This feature provides 
the computational advantage of the SES approach and 
is illustrated throughout this paper. It will be useful to 
specify an explicit one-to-one mapping between the bases 
of the associated Hilbert spaces, which we take to be 

|00---0)^ |l) = |10---0), 

|00...1)^ |2) = |01---0), 

I11---1) ^ |2«) = I00---1) . (17) 

g qubits n=2^ qubits 

The left-hand-sides are the standard computational basis 
states of the simulated g-qubit system (not to be confused 
with the computational basis states of the n-qubit SES 
processor). Similarly, any unitary quantum circuit or 
operation acting on q d-level qudits can be mapped to 
and implemented on an SES processor with n = <P qubits; 
the natural mapping is a straightforward generalization 

of (fT71) . 



matrix dimension n 


FIG. 4. (color online) Level spacing of K matrices versus n. 
Data (open circles) are averaged over 1000 random instances 
of K. The solid line is the function 1.89xn“'^’^®. The dashed 
line is the function 

The operation of a real SES chip will be nonideal, and 
it is important to consider the effects of decoherence and 
other errors on its performance. This is discussed in de¬ 
tail in Sec. IlYl The main conclusion is that although 
decoherence and unitary control errors do limit the ac¬ 
curacy of an SES computation or simulation, the effects 
of decoherence are much less restrictive here than with 
the standard gate-based approach (hence the ability to 
implement larger problem sizes). In practice, the com¬ 
plexity of fabricating a large programmable SES chip will 
likely limit its application before decoherence does. 

III. APPLICATIONS OF THE SES METHOD 
A. Uniform state preparation 

Our first example will be to generate the entangled 
state m in a single step: Consider the real nxn Hamil¬ 
tonian 

'E ~ UmaxATstarj (1^) 

where 



/I 

1 

2 

1 

2 



1 

? 

2 

0 

0 ••• 

0 

A^star — 

0 

0 ••• 

0 



0 

0 ••• 

V 


The Hamiltonian (fT31) describes a graph where qubit 1 
is symmetrically coupled to all other qubits, which are 
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where 


FIG. 5. Star network for n = 9 uniform state preparation. 


themselves uncoupled (a star network). The case of 9 
qubits is shown in Fig. [5] 

The SES chip is initially prepared in basis state |1). 
Only two eigenfunctions—let’s call them \'tjj ±)—have 
overlap with |1), so the evolution is effectively a two- 
channel problem. The spectrum is as follows: States 
have energy E± = gmax(l ± other eigenfunc¬ 

tions are degenerate with E = 0. Evolution for half a 
period corresponding to the splitting ^/ngraax^ namely 


iqc — /— 

y/ngn 

leads to the desired operation 


( 20 ) 


— il-Ltn 


11) = exp 




; ^ 
1 0 


A 

0 


Vi 0 ••• o) 


l) = |unif), 


( 21 ) 

apart from a phase. This can be implemented in a few 
ns with superconducting circuits. 

We would like to make a few remarks that apply to 
this application, as well as to many others: First, the 
magnitude of the interaction strength used in (ITOl) is ar¬ 
bitrary; any convenient interaction strength go satisfying 
0 < 5o < ffmax is sufficient. To make the operation as fast 
as possible, however, we have choosen go =ffmax- Second, 
it is not necessary to use a time-independent interaction 
strength. Any single-step “pulse sequence” of the form 


— H-Ltn 


with 


n=goK 


( 22 ) 


(23) 


and K a constant (time-independent) matrix, satisfies an 
area theorem 


jgit) dt = go tq 


(25) 


and T is the time-ordering operator. The identity (1241) 
implies that any time-dependent coupling g(t) satisfying 
(1251) can be used, simplifying experimental implementa¬ 
tion. 


B. Grover search algorithm 

Next we show how to use a programmable SES chip 
to implement the Grover search algorithm [^ . which 
introduced the powerful amplitude amplification tech¬ 
nique that has led to speedup for many other algorithms. 
Grover’s procedure for a single marked state |i) in a 
database of size n is 

(lFOi)^|unif) Ri |f), with I3=\j^\/n\. (26) 


Here 


W = 2|unif)(unif| - I 


/2 — n 

2 

2 

■■■ 2 

\ 


2 

2 — n 

2 

■■■ 2 



2 

2 

2 — n 

■■■ 2 


, (27) 

^ 2 

2 

2 

... 2- 

n/ 



is a unitary operator that performs an inversion about 
the average. 


/I 


o,= 


\ 


-1 


V 


(28) 


V 


is the oracle, a diagonal matrix with the fth element equal 
to —1 and the others equal to 1, and |unif) is the uniform 
superposition 

The W operator (071) can be implemented in a single 
step by using the SES Hamiltonian % = gmaxATfun, with 


Kfnll = 


/O 1 1 
1 0 1 
1 1 0 


1 


(29) 


Vl 1 1 ••• 0/ 


for a time 


X g-i Igit)K dt _ g- 


il-Ltn 


^qc — 


( 24 ) 


^ ^max 


( 30 ) 

















This leads to the desired operation 


|0)Hh 



/o 

I • 

• 


. TT 

1 

0 • 

• 1 


z — 

n 







I • 

■ oj 



= w, 


(31) 


up to a phase factor. 

The oracle (|28)) can be simply generated by a 27r rota¬ 
tion on qubit i. This 27r rotation can be implemented as 
a z rotation, which does not require microwaves. Each 
iteration of the amplitude amplification can therefore be 
implemented in just two steps, for any n, allowing even 
small SES chips to perform computations that would oth¬ 
erwise require thousands of elementary gates. 


C. Eigenvalue estimation 

Next we show how to use a programmable SES proces¬ 
sor to implement energy eigenvalue estimation, an appli- 
cation of the important phase estimation algorithm |65l - 
1^ that is used in many other applications. This example 
also illustrates how to translate an algorithm expressed 
in quantum circuit language to an SES protocol. 

The eigenvalue estimation procedure calculates an M- 
bit estimate of the phase (j) of the eigenvalue 
accumulated by an eigenfunction under the action of 
If the evolution time t is chosen to satisy t < 2 ttIE, 
the eigenvalue E (assumed to be positive) can be calcu¬ 
lated from E = 2TT<f>/t. To reduce the number of required 
qubits we use the iterative phase estimation circuit 
shown in Fig. IH which uses only a single ancilla. As 
the number M of desired bits of precision increases, one 
either performs a longer quantum computation—reusing 
the eigenfunction l"^)—or performs M computations in 
series, each requiring an eigenfunction preparation step. 

The algorithm measures M bits of (p one at a time, 
beginning with the least significant bit Xm, and work¬ 
ing backwards to the most significant bit Xi. Each step 
(except for the first) uses knowledge of the previously 
measured bits. We denote the bit being measured in a 
given step by m, with m = M, M — 1, M — 2, ■ ■ ■ ,1. The 
circuit for step m is shown in Fig. [6l where the rotation 
angle is 


M 

Xlm='X ^ 

J=;m + 1 

which depends on the values of the previously measured 
bits Xm-\-l 5 * * * : Xm ■ 

The main practical difficulty with prethreshold appli¬ 
cations of phase estimation is implementation of the 

controlled-^33^ 

operation, which typically requires a Trotter approxima¬ 
tion (and, in addition, a sparse Hamiltonian). However, 


|1|J) 



FIG. 6. Quantum circuit to compute the mth bit of p. Here 
H is the Hadamard gate, U = with H the model Hamil¬ 

tonian and t the evolution time, and Rz is a z-rotation. |'i/)) is 
an eigenfunction of H. The last operation is measurement of 
the first qubit in the diagonal basis; the result is Xm. G {0, 1}. 


the SES method allows any controlled unitary 

controlled—C/ (34) 

to be implemented in a single step when U can (which is 
possible when U is symmetric). To see this, assume that 
U = where A is a real NxN matrix, and write the 

2N X 2N matrix (IM)) as 

controlled—[/ = |0)(0| (E> Inxn + |I)(I|®e-*^, (35) 


where / is the identity matrix and where we take the first 
qubit to be the control. Next map the 2Wdimensional 
Hilbert space to the SES processor according to 


| 0 )®| 1 ) ^ | 1 ), 
| 0 )®| 2 ) ^ | 2 ), 


|0)®|iV)^ \N), 
|I)®|1) ^ |fV+l), 
|I)®|2) ^ |fV + 2), 


|l)G)|Af)^ |2A^). (36) 

The operation (l34l) can therefore be written as 


controlled—e 


— lA 


= exp 



(37) 


which can be implemented by an SES processor in a single 
step. The elements 0 and A on the right-hand-side of (1571) 
are each NxN matrices, with 0 the zero (null) matrix. 

We turn now to the SES eigenvalue estimation proto¬ 
col: Let H he a real NxN model Hamiltonian on which 
we wish to perform phase estimation, and denote the ba- 
sis of by {|1), |2), • • • , |A^)}. The SES implementation 
requires n = 2N qubits. The first objective in the proto¬ 
col is to prepare the initial state 


|0) ® 1^) (38) 

of Fig. [HI where |'0) is an eigenfunction of H. We 
will perform the state preparation adiabatically, which 
is restricted to states of minimum or maximum energy; 
here we prepare the ground state of El and estimate the 
ground state energy E. 
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Adiabatic ground state preparation is usually imple¬ 
mented by programming a convenient initial Hamilto¬ 
nian Hq that does not commute with H, relaxing into 
the ground state of Hq, and then slowly changing the 
system Hamiltonian from Hq to H. However, in the SES 
approach it is necessary to use nonequilibrium adiabatic 
evolution, because the physical ground state |0)®" is out¬ 
side the SES. The processor is initially prepared in the 
basis state |1). The next step is to produce the SES state 
equivalent to 


| 0 ) 0 


|l) + |2) + ... + |jV) 

Vn 


which, according to the map (1551) . is 


|l) + |2) + ... + |jV) 

Vn 


(39) 


(40) 


A possible choice when A^ is a power of two is the “trans¬ 
verse field” Hamiltonian 


log N 

Hq = —(/max ^ ^ 0"^ . (45) 

i=l 

However, the explicit matrix forms of (14511 for large N are 
complicated and the tensor-product structure is some¬ 
what artifical for our purposes. Instead we use 

^0 ~ Umax Alfuii (46) 

where ATfuu is an x matrix of the form (I29|) . The 
initial Hamiltonian (|46ll has eigenvalues 

N-l 

= -5max C'", fee {0,1,--- ,iV-l}, (47) 

i=i 


Note that (HiH) is a uniform superposition of the first half 
of SES basis states. To prepare this we use a variation 
of ([21]), namely 


_ |1) + |2) + ••• + |A^) 

y/N 


where 


'H=g„ 



(41) 


(42) 


is a 2N X 2N block-diagonal Hamiltonian. Here ATstar is 
an N X N matrix of the form (TT^ , and 0 is the N x N 
zero matrix. The operation time in (I4ip is 


where ^ The ground state is 


1 

Vn 



Vi/ 


(48) 


with energy Eq = —{N — l)gmax. The remaining eigen¬ 
functions are degenerate with energy Ei~^q = g„iax- 
At later times 0 < t < tprep the SES Hamiltonian is 
varied as 


n{t) = 


^prep ^ 

f Hq 


1 ^ 



^prep 

l 0 

oj 

1 

^prep 

1 0 

oj 


(49) 


t 


qc — 


TT 

^/max 


(43) 


This completes the preparation of the input (HUl) to the 
adiabatic evolution stage. 

At the beginning of the adiabatic evolution stage we 
program the SES Hamiltonian to be 



(44) 


where Hq is an N x N Hamiltonian with the following 
properties: 


1. Hq is real. 


Here 


max,-,/|i7„/| 

A = - ' ' (50) 

5max 

is a positive constant that ensures that j^H can be pro¬ 
grammed into the SES processor. This stage of the pro¬ 
tocol is standard: In the long fprep adiabatic limit, the 
processor will be found at t = tprep in the desired state 
(|M)l with high probability. 

Next we implement the SES equivalent of the circuit 
given in Fig. [51 beginning with the Hadamard gate 




(51) 


2. [Ho,H]^0. 

3. The ground state of Hq is the uniform superposition 
state (HO)) . 

4. The ground state is separated from the other eigen¬ 
states by an energy gap that is a nondecreasing 
function of N. 


which we write as 


where 



( sin f cos f A 
V-cosf sinf^- 


(52) 


( 53 ) 



















10 


Then we have 


H 


I InxN = — (u 0 Inxn) Kz (m'^ 0 Inxn), 
where is the 2N x 2N diagonal matrix 


(54) 


K, = 


/I 

0 

0 

0 

0 

0 

0 

0 \ 

0 

1 

0 

0 

0 

0 

0 

0 


0 

0 


0 

0 

0 

0 

0 


0 

0 

0 

1 

0 

0 

0 

0 


0 

0 

0 

0 

-1 

0 

0 

0 


0 

0 

0 

0 

0 


0 

0 


0 

0 

0 

0 

0 

0 

-I 

0 


Vo 

0 

0 

0 

0 

0 

0 

-1/ 




/ TT 

0 

0 

0 

o 

o 





0 


0 

0 

0 0 





0 

0 

TT 

0 

0 0 


— exp 

—i 

0 

0 

0 

0 

0 0 





0 

0 

0 

0 

0 





Vo 

0 

0 

0 

o 

o 



(55) 


(56) 


This leads to 




= H 0 T 


NxN, 


(57) 


where TL is the 2N x 2N Hamiltonian TL = gmaxK, with 


K = 


and 


s\v?{^)Inxn 

-COs(5)sin(f)/ArxAr 

-cos(|)sin(|) JjvxAT 

COS^(|) I NxN 


tqc — 


9n 


(58) 


(59) 


The controlled-evolution step has been discussed above 
in (l34l) through (l37l) . Applying this result to the opera¬ 
tion (|33)l leads to 


controlled—e 


= exp 


0 

0 

0 

H 


2^-h 


• (60) 


Here 0 and H are NxN matrices, with H the model 
Hamiltonian, which we assume to be real. Now let A be 
defined as in (ISOl) . Then 


controlled—e 




(61) 


where 


n = 


0 

0 

0 

{H 


and tqc = A2’” (62) 


To perform the controlled-evolution operation, the 
Hamiltonian in (|62[) is to be programmed into the SES 
processor for a time tqc- 

Finally, we implement the z rotation 


i?z(a;) 0 I NxN, 


where 


Rz{uj) = e 




g-iw/2 

0 


0 

^iuj/2 


(63) 


(64) 


This operation can be generated by applying the 2Nx 2N 
Hamiltonian H = ^max Kz for a time tqc = uj(2g^g^^. 

The final stage of the eigenvalue estimation protocol is 
the SES equivalent of ancilla measurement (see Fig. [S]), 
resulting in the observed value Xm G {0,1}. One way 
to do this is to perform a simultaneous projective mea¬ 
surement of every qubit in the SES processor. If the 
excitation at iteration m is observed to be in qubit i, we 
conclude that 


0 if 1 < i < A^, 

1 a N + l<i<2N. 


(65) 


This result follows from the correspondence (I5S1) . The 
disadvantage of this naive measurement protocol is that 
it fully collapses the SES wave function, so the eigenfunc¬ 
tion needs to be re-prepared before the next iteration. 

A simple variation of this protocol, however, avoids the 
state re-preparation step about half the time: Here we 
simultaneously measure only the first N {= n/2) qubits. 
In this case we might observe the excitation to be in qubit 
i G (I, 2, • • •, N}, or we may not find it at all. Then we 
conclude that 

{ 0 if the excitation is observed, 

1 if the excitation is not observed. 

If Xm = 0 the measurement fully collapses the state, and 
we must re-prepare the eigenfunction \ip). But if Xm = I 
we have learned only that the excitation is in the sub¬ 
space spanned by 




(67) 


which yields no information about \^p). 

It is possible to avoid the eigenfunction re-preparation 
step altogether by using an example of the ancilla- 
assisted SES method: Here we couple an n-qubit SES 
processor to an ancilla qubit, with a degree of connectiv¬ 
ity that depends on the application. (The measurement 
application requires coupling to n/2 qubits.) Alterna¬ 
tively, we can regard one of the qubits in a fully con¬ 
nected array of n -1- I qubits as an ancilla. The essential 
point is that the device now explores the single-excitation 
and double-excitation subspaces. The only disadvantage 
of this measurement protocol is that it requires N steps 
per iteration. 

The idea is to measure the multi-qubit operator 


>N ' 




NxN, 


( 68 ) 


which projects qubits i = I, 2, • • •, into a state of def¬ 
inite parity. If at iteration m the eigenvalue of (1681) is 
observed to be —I, the single excitation is in the sub¬ 
space spanned by{|I),|2),---,|A^)} and we conclude that 
Xm =0. If the eigenvalue is -|-1, the excitation is in the 
space spanned by {|A^ -I- I), |A^ -I- 2), • • •, |2A^)} and we 
conclude Xm = 1 ■ The measurement of the operator (IS51) 
can be carried out with a single ancilla qubit using the 
circuit given in Fig. [T] 

















11 



FIG. 7. Quantum circuit to measure the parity operator (I68II . 
The first qubit is an ancilla and the others are the first N 
qubits of the 2Af-qubit SES processor. The circuit uses N 
CNOT gates. 


It is useful to discuss the nonscalability of the SES 
method in the context of the eigenvalue estimation appli¬ 
cation. Typically N is exponentially large in the number 
of particles, making classical simulation impractical. An 
ideal (error-free) quantum computer would require only 
0(log N) qubits to run the phase estimation circuit of 
Fig. ini However, the large circuit depths required for the 
controlled evolutions have limited prethreshold applica¬ 
tions to very small examples. The SES implementation 
requires 2N qubits, but can perform the controlled evo¬ 
lutions in a single step. 


D. Schrodinger equation solver for 
time-independent Hamiltonian matrices 


Next we consider the problem of wave function 
propagation by a real but otherwise arbitrary time- 
independent Hamiltonian H, 

1^) ^ e-*^‘|V'). (69) 


This application, and especially its time-dependent ex¬ 
tension discussed below, play to the strengths of the SES 
chip and suggest a useful prethreshold computational 
tool. We assume that H is a real, symmetric nxn ma¬ 
trix, and we call H the model Hamiltonian. Here t is the 
length of simulated time (for example, the duration of 
some physical process). To map this problem to an SES 
processor we first find the smallest positive constant A 
such that every matrix element of 


n = 


H — const X I 
A 


(70) 


is between — ^max and gmax- Here I is the nxn identity 
matrix. When A > 1 we are “compressing” the model 
Hamiltonian down to that of the SES chip, whereas when 
0 < A < I we are expanding it. Such a rescaling is 
required because the characteristic energy scales of the 
model and SES chip are usually different. With the SES 
processor we then perform the equivalent evolution 


^ ( 71 ) 


t 


qc 


= Xt. 


(72) 


The total time required to perform a single run of the 
quantum computation is therefore 

^qu = ^qc ^meas; (73) 

where tmeas is the qubit measurement time. For super¬ 
conducting qubits we can assume tmeas to be about 100 ns 
(^ . which includes the time needed for classical post¬ 
processing. (Note that the shortest high-fidelity read¬ 
out time demonstrated to date, including resonator ring- 
down time, is closer to 300 ns l70l[ . The faster “catch- 
disperse-release” protocol of Ref. [63 has not yet been 
demonstrated.) 

A single run of the quantum computer (with read¬ 
out) simulates a single repetition of an experiment: Ini¬ 
tialization, Schrodinger evolution, and measurement. It 
is important to emphasize that such a protocol imple¬ 
ments a weak simulation, providing a single sample from 
the distribution of possible measurement outcomes, not 
the probability distributions themselves as is normally 
computed classically. (This limitation is not specific to 
the SES method and applies to state propagation with 
an error-corrected universal quantum computer as well.) 
For some applications the distinction between weak and 
strong simulation might be minor. However in other 
cases it is necessary to estimate the occupation proba¬ 
bilities pi,p 2 , ■. ■ ,Pn accurately. We discuss the runtime 
overhead for strong SES simulation below in Sec. IHI FI 

How long does a classical simulation of take? This 
of course depends on the model Hamiltonian H (includ¬ 
ing its dimension n and spectral norm), the value of t, and 
the classical processor and simulation algorithm used. To 
assess the possibility of quantum speedup, however, it is 
sufficient to find the minimum time td required to clas¬ 
sically simulate a given run of an ideal SES processor, 
with TL a “typical” SES Hamiltonian (a real nxn ran¬ 
dom symmetric matrix with all entries between —^max 
and (7max), and fqc significantly less than the coherence 
time. For this analysis we consider the case 

tqc = 100 ns and = 50 MHz. (74) 

27r 

The total quantum computation time (j73l) in this exam¬ 
ple is therefore about 

tqu = 200 ns. (SES chip) (75) 

We have studied the classical simulation runtime td 
for this problem, comparing, on a single core 0 , three 
standard numerical algorithms: 

1. State propagation via Hamiltonian diagonalization. 
For a given 'H, the unitary matrix V of its eigenvec¬ 
tors and diagonal matrix D of its eigenvalues are 
first computed. Then we numerically compute the 
product 

where l"^) is the initial state. 


where 


(76) 















12 



FIG. 8. (col or online) Classical simulation runtime on a sin¬ 
gle core versus matrix dimension n, in seconds. Here 
the computational task is solution of the Schrodinger equa¬ 
tion with a time-independent Hamiltonian by matrix diago- 
nalization. Data (circles) were determined by averaging the 
runtimes over 1000 random instances of H. The solid line is 
the function td = 203xn^'^'* ns; the scaling becomes 0{n^) at 
larger n. The runtime for a 630 x 630 Hamiltonian is about 
200 ms. 

2. Matrix exponentiation via Fade approximation 
with scaling and squaring [^ . Here we directly 
compute exp{—i'Htqc) and then multiply by \ip). 

3. Krylov subspace projection [t^. In this case the 
product exp{—i'Htqc)\ilj) itself is directly calculated. 

In all cases we assume an initial state |'0) of the form 

0 


which corresponds to a single SES basis state, and we av¬ 
erage the computation times over 1000 random instances 
of TL. Although the three methods have similar speed 
and accuracy for the particular problem simulated here, 
the matrix diagonalization method was the fastest, fol¬ 
lowed by matrix exponentiation. We also tested Runge- 
Kutta integration and matrix exponentiation via Cheby- 
chev polynomial expansion , which were not compet¬ 
itive with the above methods for the specific application 
considered. In Fig. |S] we plot the measured single-core 
runtimes for the optimal classical algorithm (Hamilto¬ 
nian diagonalization) versus matrix dimension n. We 
observe that the quantum simulation time (1751) is much 
shorter than all of the single-core runtimes considered. 


Our objective is to achieve speedup relative to a state- 
of-the-art supercomputer, not a single core. The classi¬ 
cal simulation runtime tci should then be evaluated on 
a supercomputer, using an optimally distributed parallel 
algorithm. However, we can bound the parallel perfor¬ 
mance by using the single-core result and assuming per¬ 
fect parallelization efficiency: We approximate a petaflop 
supercomputer by 10® gigaflop cores, and conclude that 
the classical runtime can be no shorter than 10“® times 
the single-core time. (This is a conservative estimate be¬ 
cause high parallelization efficiency is not expected for 
problem sizes smaller than the number of cores.) We 
therefore conclude that, for this particular state propa¬ 
gation application, the classical simulation runtime is no 
shorter than 

tci = 203xn^'^‘*fs, (classical supercomputer) (78) 

while the quantum simulation can be performed in a few 
hundred nanoseconds. The breakeven dimension accord¬ 
ing to d78l) is about n = 630 qubits. This is quite large 
given the full connectivity requirement, and it is not 
known whether such a device could be built in practice. 
However the breakeven dimension in the time-dependent 
case (discussed below) is considerably smaller. 

In our estimate of the classical simulation runtime we 
have not included the time needed to store the Hamil¬ 
tonian matrix in memory or perhaps compute it from a 
separate procedure. Similarly, for the quantum simula¬ 
tion time estimate we have not included the time required 
to send the n(n-|-l)/2 controls to the qubits and couplers 
before the simulation. Furthermore, not every n > 630 
simulation will exhibit a speedup; this depends on the 
particular simulated Hamiltonian and the simulated time 
duration t. 

An interesting aspect of the Schrodinger equation 
solver is that the complexity is 0(1): The quantum sim¬ 
ulation time is independent of n. This implies that the 
SES method yields an exponential speedup for this ap¬ 
plication. However such complexity considerations are 
probably not meaningful given that the method is not 
scalable. 


E. Schrodinger equation solver for time-dependent 

Hamiltonians: Simulation of molecular collisions 

Finally, we discuss what is perhaps the most interesting 
application of the SES method known to date, the solu¬ 
tion of the Schrodinger equation with a time-dependent 
Hamiltonian. This is a straightforward generalization 
of the time-independent case, but we expect the time- 
dependent case to be more useful in practice. In this 
section we provide a detailed example of time-dependent 
Hamiltonian simulation with a small SES chip. 

Time-dependent Hamiltonian simulation is imple¬ 
mented by varying the SES matrix elements © accord¬ 
ing to some protocol, which can be done with nanosecond 
resolution. This does not require any additional runtime. 
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the time complexity is still constant, and the total quan¬ 
tum simulation runtime for a 100 ns evolution is again 
given by dZSl). Although the classical runtime is problem 
specific, we can again assess the possibility of speedup 
by estimating the time required to classically simulate 
an ideal SES processor, in this case with all matrix el¬ 
ements varying on a nanosecond timescale. There are two 
types of numerical simulation algorithms we consider: 

1. Runge-Kutta integration. Here we solve the system 
of coupled ordinary differential equations 

d = —iHa. (79) 

Although the Runge-Kutta runtime is slower than 
diagonalization for a time-independent Hamilto¬ 
nian, it does not slow down significantly when T-L 
is time dependent. 

2. Time slicing combined with diagonalization. This 
algorithm is based on an approximate decomposi¬ 
tion of the time-dependent problem into a sequence 
of constant-Hamiltonian intervals, each of width 
At. The time At must be significantly smaller than 
the characteristic timescale of matrix element vari¬ 
ation, for example At = 0.1ns. Then 

iV^lice = ^ (80) 

time slices are required, and the classical runtime 
using this approach will be approximately JVsHce 
times longer than (1751). For the 100 ns evolution, 
A^,lice=1000. 

We find that Runge-Kutta integration is the fastest ap¬ 
proach for the specific problem considered here. In Fig.^ 
we plot the measured single-core runtimes for the opti¬ 
mal classical algorithm (Runge-Kutta integration) versus 
matrix dimension n. Bounding the performance of this 
algorithm on a petaflop supercomputer by including a 
factor of 10“® (recall discussion from Sec. lHIDl . we con¬ 
clude that for this application the classical simulation 
runtime is no shorter than 

tci = 1.38xn^'^®ns, (classical supercomputer) (81) 

while the quantum simulation can be performed in a few 
hundred nanoseconds. The breakeven dimension accord¬ 
ing to is around n = 50 qubits. As expected, this 
value is much smaller than the breakeven in the time- 
independent case. An SES Schrodinger equation solver 
of modest size might be able to achieve quantum speedup 
relative to a petaflop supercomputer. 

We turn now to a detailed example of time-dependent 
Hamiltonian simulation with a small programmable SES 
chip. One particularly interesting application of the 
method is to the quantum simulation of atomic and 
molecular collisions. Collisions are especially well suited 
for SES simulation because they typically involve mod¬ 
est Hilbert spaces—tens to thousands of channels—and 



FIG. 9. (col or online) Classical simulation runtime on a sin¬ 
gle core versus matrix dimension n, in seconds. Here 
the computational task is solution of the Schrodinger equa¬ 
tion with a time-dependent Hamiltonian by Runge-Kutta in¬ 
tegration. Data (circles) were determined by averaging the 
runtimes over 1000 random instances of R. The solid line is 
the function td = 1.38 x ms. The simulation time for a 
50 X 50 Hamiltonian is about 200 ms. 


in the time-dependent formulation involve Hamiltonians 
that are naturally bounded in time. In particular, the ini¬ 
tial and final asymptotic Hamiltonians for neutral scat- 
terers are diagonal (in the adiabatic basis), whereas the 
off-diagonal elements rapidly turn on and then off dur¬ 
ing the collision itself, inducing transitions between the 
channels. Although the Born-Oppenheimer potential en¬ 
ergy surfaces used here do require a classically inefficient 
electronic structure precomputation, the largest poten¬ 
tial energy surface calculations [zl, [t^ are far ahead 
of the largest classical collision simulations performed 
to date [TTI - IT^ . SES implementation of the semiclassi- 
cal Born-Oppenheimer problem therefore has the poten¬ 
tial to push molecular collision simulations to new unex¬ 
plored regimes. (We note that there are related chem¬ 
ical reaction simulation methods developed by Lidar et 
al. and by Kassal et al. [8l| that do not require a pre¬ 
computed potential surface, but these require an error- 
corrected quantum computer to implement and are not 
prethreshold methods.) Another useful feature of the 
scattering application is the convenient mapping of each 
molecular channel to a single SES basis state, which is 
possible because of the similar way initial states are pre¬ 
pared and final states measured in both the processor 
and a collision experiment. We also find that the atomic 
physics time and energy scales turn out to map nicely to 
that of superconducting qubits after optimal rescaling. 

To illustrate this application we consider a three- 
channel Na-He collision (an unpublished preliminary ac- 
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FIG. 10. (color online) Time dependence of the matrix el¬ 
ements of the scattering Hamiltonian for a collision with 
Vo = 2.0 and b — 0.5 in atomic units. The diagonal matrix 
elements (solid curves) are similar in magnitude and cannot 
be resolved in this figure. The off-diagonal elements (dashed 
curves) are much smaller than the diagonal elements and also 
cannot be resolved here. In this example the collision energy 
is Ecm ~ 341 keV and the impact occurs at to = 6 x 10“^® s. 

count of this application is given in Ref. [s^ ). The three 
channels included in our model and their correspondence 
with SES basis states are 

Na(3s)-bHe(ls2) [1 2S+] ^^ |l) (82) 

and 

Na(3p)-bHe(ls2) [1 2n;2 2 e+] ^ |2),|3). (83) 

The square brackets indicate the molecular structure of 
the channels. In this model, the helium atom remains in 
its electronic ground state Is^ during the collision (the 
excitation energies of its excited states are too high to 
be relevant here), whereas sodium can be excited from 
its ground state 3s to either of two excited states, both 
denoted by 3p. In the physical system, the channels (15^ 
and (l83l) have additional degeneracies, including spin de¬ 
generacies, but they do not affect the collision probabil¬ 
ities calculated here. Precomputed Born-Oppenheimer 
energies and nonadiabatic couplings of the Na-He system 
are stored for fixed values of the internuclear distance 
R, and we make a standard semiclassical (high energy) 
approximation and assume that the scatterers follow a 
straight-line trajectory. Then, for an impact occuring at 
a time the internuclear separation varies according to 

R = IP + vl{t - toY, (84) 

where Uq is the initial relative velocity and h the impact 
parameter of the collision. The relative velocity is re¬ 
lated to the collision energy in the center-of-mass frame 



FIG. 11. (color) Scattering probabilities pi_>i(t) for a Na- 
He collision with vo = 2.0 and b — 0.5 in atomic units. The 
system is initially prepared in channel 1. The collision occurs 
at to = 6x 10“^® s. 

through Ecm = Mfo/2, where /i is the reduced mass. The 
procedure outlined in Appendix ICl then leads to the scat¬ 
tering Hamiltonian shown in Fig. HOI 

In Fig.[TT]we plot the probabilities pi^i(t) for the Na- 
He system to be found in channel i S {1, 2, 3} after being 
initially prepared in channel 1, the ground state. The 
final values pi_>i(oo) are the probabilities for an elastic 
(i = 1) or inelastic (i = 2,3) collision with a given Ecm 
and 5, which we find to be 

Pi^i = 0.116 

pi-^2 = 0.038 

pi^3 = 0.846. (85) 

To simulate this process with a programmable SES 
chip we must hrst rescale the physical or model Hamil¬ 
tonian so that it hts on the SES processor. Doing this 
optimally is critical to the utility of any time-dependent 
simulation so we will discuss it in some detail: Suppose 
for the moment that our model Hamiltonian is given by 
a time-independent nxn real symmetric matrix H. Di¬ 
viding E[ by any positive constant A while rescaling the 
evolution time by the same factor obviously leaves the 
dynamics invariant. Because we want the quantum sim¬ 
ulation to be as fast as possible, we choose the smallest 
value of A that makes H/X compatible with the SES pro¬ 
cessor (every matrix element of E[/X is between —(/max 
and (/max)- As mentioned above, if A > 1 we are com¬ 
pressing the model’s energy scales to fit on the SES chip, 
whereas if A < 1 we are expanding them. This naive ap¬ 
proach to rescaling, however, does not take advantage of 
the fact that we can always shift H by a constant (which 
changes the corresponding states by a phase factor that 
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FIG. 12. (color online) Rescaling function for the Na-He colli¬ 
sion simulation with (;inax/27r = 30 MHz. Collision parameters 
are vo = 2.0 and 6 = 0.5 in atomic units. We find that A has an 
asymptotic value around 10^ and reaches 6 x 10® during the 
collision. 


FIG. 13. (color online) Nonlinear time scaling for the A(t) 
given in Fig. 1121 Most of the simulation time is spent near 
the moment of collision, to, and a single run of the simulation 
is completed in 24 ns. 


we do not measure). Including this gauge transformation 
results in the rescaling used above in (1701) . The time tqc 
required to simulate an evolution of duration t is sim¬ 
ply given by ([721), but now we will go further and regard 
(I72p as giving the linear relationship between the simu¬ 
lated and physical times during a process. To generalize 
this construction to a time-dependent model Hamiltonian 
H(t), we first compute the mean of its diagonal elements, 

1 ” 

c{t) = -y^Hu, ( 86 ) 


and then hnd, at each time t, the smallest positive A such 
that every matrix element of 


m) 


H{t) — c{t) X I 

W) 


(87) 


is between —^max and (/max- The function X{t) defines 
the resulting nonlinear relation between the physical and 
simulated times according to 

tqc{t) = [ ( 88 ) 

Jo 

Equation (1871) gives the simulated Hamiltonian as a func¬ 
tion of the physical time t, and ((551) is then inverted to 
find the desired ’H(tqc), which in turn is programmed into 
the SES chip. 

Applying this procedure to the Na-He collision prob¬ 
lem results in the rescaling function A shown in Fig. 1121 
and the nonlinear time relationship shown in Fig. 1131 We 


note that the nonlinear energy/time rescaling protocol is 
extremely effective at mapping this atomic physics prob¬ 
lem to the SES processor, allowing a single run of the 
simulation (excluding measurement) to be completed in 
about 24 ns. The A function shown in Fig. [T51 assumes 
ffmax/27r = 30 MHz; if this is increased to 50 MHz the 
simulation is completed in 15 ns. It is important to em¬ 
phasize that any positive piecewise continuous function 
A(t) defines a mathematically valid energy/time rescal¬ 
ing, and that the specific form used in practice should 
be determined by hardware considerations, such as qubit 
coherence times and control pulse bandwidth. In partic¬ 
ular, X(t) can be chosen to bound both the magnitude of 
the SES matrix elements and their rates of change, but 
we will not pursue this variation here. 

Use of the rescaling function given in Fig. ITT] leads to 
the SES matrix elements shown in Fig. 1141 which bear no 
resemblance to those of Fig. [TUI The corresponding scat¬ 
tering probabilities during the simulation are shown in 
Fig. [151 Compared with Fig. [TTl we see that the dynam¬ 
ics near the moment of collision are rescaled to occupy 
most of the simulation. The final scattering probabilities 
are the same as in Fig. ITT] and are given in ((55|) . 


F. Strong quantum simulation 

As we have emphasized, a single run of the SES chip 
provides a single sample from the distribution of possi¬ 
ble measurement outcomes—a weak simulation—not the 
probability distributions themselves as is normally com¬ 
puted classically. If the objective is to perform a strong 
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FIG. 14. (color) SES Hamiltonian matrix elements for a Na- 
He collision with vo = 2.0 and b — 0.5, in atomic units. Here 
we assume that yniax/27r = 30MHz. At each instant the mag¬ 
nitude of at least one matrix element achieves its maximum 
value of 30 MHz, making the simulation as fast as possible. 



FIG. 15. (color) Scattering probabilities in the SES 

processor for a Na-He collision with vo = 2.0 and b — 0.5 in 
atomic units. The dynamics near the moment of collision 
occupy most of the simulation. 


simulation and hence estimate the n basis state occupa¬ 
tion probabilities pi , it is necessary to repeat the simula¬ 
tion and readout many times. In this section we dis¬ 
cuss the runtime overhead for strong SES simulation. 
(The objective considered is that of measuring occupa¬ 
tion probabilities, not probability amplitudes.) 

Suppose that after some simulation we want to mea¬ 


sure the occupation probability p of one qubit in the SES 
chip. We do this by performing the simulation N times, 
after each repetition r measuring the qubit in the diago¬ 
nal basis and observing x S {0,1}. The estimate 

1 ^ 

Pest = ^ ^ (^9) 

r—1 

for p will have a sampling error (standard error of the 
mean) given by 

E ^ ^var(pest) = y(90) 

Eor example, to ensure that the sampling error is smaller 
than 1%, it is sufficient to repeat the simulation 2500 
times. 

The overhead N = 2500 is a worst-case estimate for a 
1% sampling error. If p is known to be small (or close to 
1), fewer repetitions are required. However in this work 
we restrict ourselves to the case where measuring p is the 
objective of a quantum simulation and is not known a 
priori. 

We turn now to the runtime overhead for estimat¬ 
ing every SES basis state occupation probability pi in 
an n-qubit processor. The qubits are measured simul¬ 
taneously and the sampling error formula (j90|) apples to 
each Pi. The only modification resulting from the SES 
constraint—that the device is in the single-excitation 
subspace—is that the n probabilities are not indepen¬ 
dent, because 


= (91) 

i=l 

However the condition (l?T]) does not affect the sampling 
statistics. Therefore we conclude that the sampling error 
for the qubit is 


E,; = 


N 


(92) 


The result (IMll shows that the strong simulation over¬ 
head, given by the required number of repetitions N, is 
independent of n and is no worse than that for a single 
qubit. In particular, the upper bound 
applies, and hence the required number of repetitions 
satisfies 


N < 


4E? 


(93) 


In conclusion, the complexity for strong quantum simu¬ 
lation is also constant^ and the overhead factor N is no 
larger than the single-qubit value. 

Although the sampling error can be made arbitrarily 
small, it is usually not helpful to require it to be smaller 
than the deviations in the occupation probabilities re¬ 
sulting from other error sources, such as decoherence. 
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Consider, for example, the evolution of a single qubit for 
100ns in the presence of Ti relaxation. If Ti =40p,s, deco¬ 
herence would lead to a 0.25% error in the excited state 
probability p. It is not possible, in general, to measure 
the original qubit excitation probability with better than 
99.75% accuracy because of this error. Thus, decoherence 
limits the accuracy of a strong simulation independently 
of the sampling error. We will discuss decoherence and 
other error sources in detail in the next section. 


and 

E = (100) 

Ti 

which is independent of n. An SES state (in the absence 
of unitary evolution) thefore relaxes at the same rate as a 
single excited qubit. The approximation in (llOOl) applies 
when tqc Ti, which is the regime of interest here. 


IV. ACCURACY OF SES COMPUTATION 


B. Pure dephasing error 


In this section we discuss the errors incurred during a 
quantum computation or simulation with an n-qubit SES 
chip. We are especially interested in the n-dependence of 
these errors, for large n, and whether they pose any seri¬ 
ous limitation to the practical utility of the SES approach 
(we conclude that they do not). Below we separately ana¬ 
lyze decoherence errors, matrix element “control” errors, 
and leakage out of the SES. In each case the ideal, error- 
free state after some process is an SES pure state IV'ideai), 
and we estimate the error 

E = 1 ('^ideal IpIV^ideal) 5 (^^) 

where p is the realized final state. 


Next we discuss pure dephasing, which in Xmon qubits 
is believed to be caused primarily by flux noise (how¬ 
ever this has been recently questioned (H)- We again 
estimate this error by setting 'Hw — 0 and calculating 
the degradation of an initially prepared SES state 
in the absence of unitary evolution. We assume a stan¬ 
dard single-qubit dephasing model with no correlations 
between the noise at different qubits. 

The Kraus matrices in this case are 


El 



and E 2 



( 101 ) 


with 


A. Energy relaxation error 


— d 2iqc/T,/, 


( 102 ) 


The first decoherence error we discuss is energy relax¬ 
ation (zero-temperature amplitude damping). We esti¬ 
mate this error by setting T-Ln' = 0 and calculating the 
decay of an initially prepared SES state 

n n 

l^/’init) = ^ a* |i) = ^ ttj |0 • • • li • • • 0) (95) 

i—1 i—1 

in the absence of unitary evolution. 

The single-qubit Kraus matrices for this process are 

= (0 X) “<> a = (S ■ (“O) 

with 

r = e"*‘>=/^i. (97) 

Here tqc is the runtime for some SES quantum compu¬ 
tation and all qubits are assumed to have the same Ti 
value. The Ti time for capacitively coupled Xmon qubits 
is currently about 40 ps . 

In the presence of energy relaxation, 

Pinit = |V^i„it)(V'init| ^ 4®” + • • • , (98) 

where the dots denote terms involving one or more ap¬ 
plications of the E 2 operator ( 1 ^ that are outside of the 
SES and do not contribute to (IMl) . Then 

(99) 


We can estimate the time for Xmon qubits (with fixed 
capacitive coupling) by using the relation 


1 _ 1 1 

%~Y2~Wi^ 


(103) 


with values Ti = AO ps and T 2 — 20 ps from Ref. [^ . 
which leads to Ri 27 ps. (Note, however, that the flux 
noise in Xmon qubits is not Markovian, as we have as¬ 
sumed. We believe that our simple dephasing calculation 
overestimates the actual dephasing error.) All qubits are 
assumed to have the same value. 

In the presence of pure dephasing, 

Pinit rLi Pinit r^i 
n 

-E ^ (El® • • ■ ■ ®Ei) Pinit (El® • • • ^2 • • • ®Ei)\ 

qubit i 

where the Ei ^2 now refer to (IIOII) . and we have used 
the fact that terms with two or more applications of the 
E 2 operator vanish when applied to pinit- Note that E 2 
annihilates an SES basis state |i) = |0 • • • li • • • 0) unless 
it acts on qubit i, in which case it produces a factor of 
i/l — r. Then dephasing transforms an SES state to 


n 

Pinit r Pinit -I- (1 - r)^ |oip |i)(j|, 
2=1 


Pinit ^ Pinit + ' ' ' 


( 104 ) 
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and the associated fidelity loss is 
E = (1 - 

(105) 

^ i=i ^ 

The dephasing error (jl05l) is maximized when the SES 
basis states are equally populated, |ai| = l/-v/n. Then 

max E = (1 — 6“^*'’°/^'^) (l -^ ~ ^ (1 -^ . (106) 

at \ n J \ n J 

This expression is valid for n > 1 (the n = l SES state 
|1) has no pure dephasing error). The n-dependence of 
the worst-case dephasing error (11061) is very favorable, 
approaching a constant as n —> oo. Therefore the pure 
dephasing error is bounded by 

E < 1 - (107) 

which is only a few times larger than (HOPE 

The total error due to decoherence is the sum of (jlOOIl 
and (|105l) [or (11071) ]. Assuming a tqc = 100 ns computa¬ 
tion and the coherence times given above, this error is 
around 1% and is independent of n. 

C. Hamiltonian control errors 

Next we calculate the error (iMl) caused by imper¬ 
fect experimental programming of the SES matrix ele¬ 
ments, which we call a control error. We assume that 
the intended Hamiltonian is a real, symmetric, time- 
independent nxn matrix 77, but that the applied Hamil¬ 
tonian is instead % + V, where H is a real, symmetric, 
time-independent matrix that does not commute with 'H . 

We consider a typical situation where the processor is 
initially prepared in a single SES basis state |i). The 


nonideal final state (neglecting decoherence) is then 

(108) 

where t is the evolution time for either a complete algo¬ 
rithm or a single step in an algorithm. The error (1941) in 
this case is therefore 

E, = 1 - (109) 

Averaging (jlQ9p over the initial SES basis state leads to 

1 ^ 

E = 1 - - y |(7|e*^*e-*(^+^)‘|7)|^ (110) 

n ' ' 

i=l 

The SES Hamiltonian 77 in (jllOl) is assumed to have the 
“typical” form described above in Sec. Ill Cl 
In this section we evaluate the control error (jllOl) using 
two complementary approaches. First we consider the 
small-H7 perturbative limit. Evaluating quantities of the 
form 

by a series expansion in H is a standard problem in per¬ 
turbation theory: Differentiating (Hill) with respect to 
time yields 

r)JT 

i— = V{t) U, with V{t) = (112) 

showing that Hlip satisfies a Schrddinger equation with 
time-dependent Hamiltonian V. We can therefore write 

(mm as 

u = (113) 

where T is the time-ordering operator. Expanding H13I) 
to second order yields 


U=l-if dh V{ti) - \ fdt, fdt2 T{V{h)V{t2)) 
Jo ^Jo Jo 


= l-i [ dtiV{ti)-l-[ dh{ I dt2Vih)Vit2)+ I dt2V{t2)Vih)]. 


(114) 


This leads to 


{i\U\i) = l-i I dti - ^ / dh / dfa |i) 


+ [ =1-*E f dhS^J(-h)Sk^{h)VJk 

Jti J -i. Jo 

- ^ E / d7i f dt2S^j{-h) Skiih-t2) Smiih) + j dt2 Sij{-t2) Ski(t2 - h) Smi{h) (115) 

^ jklm L7o Jti 


































19 


where 


Sii'{t) = {i\e 


(116) 


Because % is symmetric, Su'{t) is also a symmetric matrix. Then to second order in V we have 

|(i|C/|z)|^ = 1 + 2Im^ / dtiSij{-ti)Ski{ti)Vjk+'^ [ dti f dt[ Sij{-ti) Skiiti) Sl,^{t[)VjkVj 

,1. JO j^kk' “^0 -^0 


jk 

nt 


- Re ^ f dti f dt 2 Sij{-ti)Skl{tl-t 2 )Smiit 2 )+ f dt 2 Sij{-t 2 )Skl{t 2 -ti)Smiiti) 
Ldo dti 

The control error to second order for a fixed H and V follows from (IllOl) and (11171) . 

I- 


VjkVlm. (117) 


Next we regard the perturbation as random and aver¬ 
age over the random matrix V. The elements Vi<ii are 
assumed to be independent identically distributed ran¬ 
dom variables (the elements Vi>i' fixed by symmetry). 
The moments of V are 

'Wb = 0 (118) 


and 


example is shown in Fig. 1161 for t = 10ns and matrix ele¬ 
ments of the perturbation V uniformly distributed in the 
interval 


J’ 


( 122 ) 


with 5VI2'it = 0.5 MHz. The value of tr for this distribu¬ 
tion is 


^ah^cd — ^ {dacdbd T dad^bc dac^bd^ab^ ■> (H^) 


where cr is a parameter with dimensions of energy char¬ 
acterizing the size of the control errors. The condition 
(I119|) enforces the symmetry requirement of V. The 
perturbation-averaged control error (IllOl) is then found 
to be 


9 9 /"* 

E = -\ -/ dti dt2 


TvS{ti-t2)\ 

n 

- 2^ \Su{tl - t2)|' - 51 + ^2)1 


2=1 


+ E i^b(ii)r i^b(i2) 


2=1 

|2 


( 120 ) 


where the propagator Su' is defined in (|116l) . The con¬ 
trol error is proportional to cr^, as expected, and for fixed 
evolution time t and large n is dominated by the second 
term (the first term in the square brackets). By retaining 
this dominant second term, together with the first 
term, and evaluating the trace in the eigenfunction ba¬ 
sis, we can perform the time integrations analytically to 
obtain a useful spectral form for the perturbative control 
error, 


E 


2ah^ + — y 

n ^ 


a, OL 
Ol^Ol 


1 — COs[{Ea — Ea')t] 

{Ea.-Eadf 


( 121 ) 


where the Ea are the eigenvalues of %. 

The expression (11211) is useful for studying the n- 
dependence of the control error for short times, corre¬ 
sponding to a single step in an SES computation. An 


dV 

~ Vl2 ■ 


(123) 


We conclude from Fig. [16] that the expression (I121|) pro¬ 
vides an accurate approximation for the short-time con¬ 
trol errors, and that they are less than 0.1% in all the 
cases considered. 

The perturbative result (11201) can also be used to un¬ 
derstand the long-time asymptotic limit of the control 
error. For fixed dimension n and large evolution time t, 
we find that 


E Ri 2a^f, (124) 

independent of n. Both the first term in (11201) and 
the second term contribute to (|124l) . From (11241) we ob¬ 
tain an upper bound 


t 


max 


1 

V2a 


(125) 


on the evolution time t for which perturbation theory is 
valid. 

It is also useful to directly calculate the control error 
(I110|) numerically, which is useful for exploring the non- 
perturbative regime. The results are shown in Figs. 1171 
and m The most important conclusions of these simu¬ 
lations is that the dimension dependence of the control 
error grows much more slowly than linearly, and for errors 
in the SES matrix elements < 0.25 MHz is about 2% when 
n = 100 and 3% when n = 1000. If the size of the errors 
in the matrix elements are doubled—to <0.50 MHz—the 
n = 100 control error increases to 7.4%, less than that 
predicted by the quadratic scaling resulting from per¬ 
turbation theory. 
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matrix dimension n 



matrix dimension n 


FIG. 16. (color online) Plot of the control error E versns ma¬ 
trix dimension for gmax/Syr = 50 MHz, t = 10 ns, and errors in 
the individual matrix elements uniformly distributed between 
-0.25 and 0.25 MHz. Data (open circles) follow from an exact 
numerical calculation of (ITTOt . averaged over ideal Hamilto¬ 
nian T-i and perturbation V. The solid curve is the function 
E = 8.1 X 10“® X The dashed curve follows from the 

expression (HSJ. 



matrix dimension n 


FIG. 17. (color online) Fidelity loss caused by control errors 
in the SES Hamiltonian versus number of qubits or matrix 
dimension n. Here we assume (?max/27r = 50 MHz, t=100ns, 
and errors in the individual matrix elements uniformly dis¬ 
tributed between -0.25 and 0.25 MHz. Errors (open circles) 
are averaged over ideal Hamiltonian T-i and perturbation V. 
The solid line is the function E = 8.2x 10“® 


FIG. 18. (color online) Same as Fig.[T7]for larger SES matrix 
dimension. Here (7max/27r = 50 MHz, t=100ns, and errors in 
the individual matrix elements uniformly distributed between 
-0.25 and 0.25 MHz. The solid line is the function E = 5.7 x 


D. Leakage out of the SES 

We briefly comment on two additional error mecha¬ 
nisms that result in excitation out of the SES, namely 
leakage into the triple-excitation subspace and |2)-state 
errors. To understand the origin of leakage into the 
triple-excitation subspace, we use the identity [see ([3])] 
= c -b c’f and expand the qubit-qubit interaction in 
Terms proportional to c|ci' and Cicj, connect SES 
states to other SES states, whereas the terms propor¬ 
tional to cjcj, connect SES states to triply excited states. 
However these excitations are protected by large energy 
gaps and the corresponding errors are negligible. 

To understand the origin of |2)-state errors, recall that 
with conventional gate-based superconducting quantum 
computation, the dominant sources of |2)-state excitation 
are microwave pulses and two-qubit gates. However nei¬ 
ther of these are used in the SES approach. Initialization 
of SES basis states can produce |2)-state errors, but it is 
known how to limit these errors to less than 10“^ 1^. 


V. CONCLUDING REMARKS 

The SES method described here appears to be dis¬ 
tinct from previously investigated approaches to quan¬ 
tum computation. Like analog quantum simulation, 
it is tied to a specific hardware model and cannot be 
implemented on any architecture. However, the ap¬ 
proach enables universal quantum computation and sim¬ 
ulation, and might make quantum speedup possible with 
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prethreshold hardware. 
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Appendix A: General qubit-qubit coupling types 

In this section we discuss the generalization of the SES 
method to fully connected quantum computer models of 
the form 

Hqc = 'y ^ 4- — y ( gw y ( ® cTj/, (Al) 

i ii' 

where the (with fj, = x,y, z) are Pauli matrices and 
is a fixed, real, dimensionless tensor. In this case the 
SES matrix elements are 


\Hqc\i') = 

^ 9ij ) Jzz + ( 9jj' ) Jzz 

^ 3 3 < 3 ' 

+ 

Jxx Jyy '^(^Jxy Jyx^ 

gii'- 


It is interesting to notice how decoherence only barely 
limits the problem sizes that can be implemented with 
a programmable SES processor, at least for the applica¬ 
tions explored here. The main factor limiting the utility 
of the SES method is the difficulty of building fully con¬ 
nected arrays of qubits. In this sense we can say that 
the SES approach trades the familiar limitations result¬ 
ing from decoherence for a new limitation—that of build¬ 
ing qubit graphs with high connectivity. Relative to the 
large, community-wide effort devoted to studying and im¬ 
proving quantum coherence, the problem of increasing 
connectivity is certainly in its infancy. 

If we accept that the SES method outperforms the tra¬ 
ditional gate-based approach for prethreshold universal 
quantum computation, but that it is ultimately unscal¬ 
able, the question becomes whether an SES chip could 
be built that is large enough to be of practical use (and 
before an error-corrected universal quantum computer 
arrives). We speculate that for selected applications, 
breakeven with a classical million-core supercomputer is 
possible, a significant feat, but that a stronger form of 
speedup—such as performing computations that are es¬ 
sentially impossible classically—is probably not. Perhaps 
a breakeven-sized SES processor would be useful for its 
low power consumption or as a special-purpose subunit 
in a conventional error-corrected quantum computer. 


Note that the term proportional to 9jj’)^ii' is an 

energy shift and can be dropped. The SES method can 
be applied (possibly with some protocol modifications) 
whenever 

Jxx + Jyy ^ 0 (-^ 3 ) 

and 

Jxy — Jyx • ('^^) 

The condition (jAjp means that the interaction has an 
exchange or transverse component, and (IA4I) ensures that 
the SES Hamiltonian is purely real. 

Appendix B: Tunable coupler circuit 

In this section we calculate the qubit-qubit interaction 
strength g for the coupler circuit shown in Fig. 1191 which 
is the building block for a programmable SES chip. We 
first give a simplihed treatment by making weak coupling 
and harmonic approximations, and then discuss the gen¬ 
eral case afterward. 

The circuit of Fig. [19] has six active nodes (black dots) 
and is described by six node flux coordinates [s^. How¬ 
ever, all nodes except those labeled Ui _2 have negligi¬ 
ble capacitance to ground and are therefore “massless” 
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Vj V2 



FIG. 19. (Color online) Basic conpler circuit for SES proces¬ 
sor. The three crosses represent Josephson junctions, each 
with a fiux-tnnable Josephson coupling energy. Each Xmon 
qnbit has capacitance C and tnnable junction indnctance Lj. 
The coupler wire has a Josephson junction with tunable in¬ 
ductance Lc and mutual inductance m as indicated. Self in¬ 
ductances of the coils are denoted by Lq and Lq. 



FIG. 20. (Color online) Replacing the inductive network of 
Fig.Hniby an effective circuit with mutual inductance M and 
self inductances Lq. 


degrees of freedom that remain in their instantaneous 
ground states. They will be eliminated from the problem 
in the analysis below. The Lagrangian for the circuit of 
Fig. [His 



where $o = ^/2e is the flux quantum, C is the qubit ca¬ 
pacitance, is the dimensionless node flux at Vi, 2 , and 
U is the total potential energy. Following the approach 
of Ref. [6l|, we replace the inductive network of Fig. [T9| 
(excluding the capacitors) by the equivalent circuit of 
Fig. IH where M and Lq are effective inductances to be 
determined in terms of the physical circuit parameters. 

The potential energy in the circuit of Fig. [His Lq(/J-|- 
/|)/2 -I- MJ 1 / 2 . Rewriting this in terms of magnetic flux 


we have 


^ /Lq M\ flA 
{^ 2 ) \M lJ {I 2 ) 


u = 


(J )2 

2KL„ 


2KL, 


-I- Fii $1$2, 


(B2) 


(B3) 


where 


f M\'^ M 

The cross term in (IB3I) proportional to Fn is responsible 
for the qubit-qubit coupling. In the weak coupling limit, 
M <C Lq and LT Ri 1, which we assume below. 

Next we calculate the qubit-qubit interaction strength 
g induced by this cross term. Perhaps the simplest way 
to do this is to use the expression 




for the flux of an LC oscillator with frequency 

I 


(B5) 


(B6) 


in terms of creation and annihilation operators. This 
leads to the desired result [6l| 


9 = 


FllLq 

2 


e. 


(B7) 


Finally, we find M and Lq in terms of the physical 
circuit parameters. Assuming an time dependence 
we have from (IB2I) that 


1 


M = —X , 

iijj V ^2 


Ri 


and Lq = 


//i=0 

Using these expressions we find 


iio V Ii 


(B8) 


Jo—0 


M = 


2L 


j and Lq = Lj -I- Lq — M, (B9) 


from which we obtain 


g = - ^7(BIO) 

2(Lj J- Lo)(Lc J- 2Lg) 

Here e is the qubit frequency (jB6l) . This expression for 
the strength of the transverse coupling in the 

weak coupling and harmonic approximations is the main 
result of this section. The advantage of the circuit of 
Fig. IH over that of Ref. when extended to many 
qubits is the absence of coupler loops through which the 
flux must be individually controlled. 

Based on our recent work [? ] on a closely related cou¬ 
pler circuit, we expect the result (IB10|) to be a good ap¬ 
proximation to the actual coupling. The main difference 
is that the qubit anharmonicity suppresses the magnitude 
of the coupling (in Ref. [? ] the coupling was found to 
be suppressed by about 15%). Anharmonicity also gener¬ 
ates a small (< 1 MHz) diagonal interaction, but 

such an interaction has no effect on a single excitation. 
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R (m) x10'''° 


FIG. 21. (color online) Diagonal elements Ha of the Na-He 
scattering Hamiltonian with collision parameters no = 2.0 a.u. 
and 6 = 0.5 a.u. The three curves cannot be resolved in this 
figure. 

Appendix C: Scattering Hamiltonian 

In this section we outline the construction of the scat¬ 
tering Hamiltonian for the Na-He collision discussed in 
Sec. IIII El We begin by constructing a 3 x 3 matrix U as 
a function of internuclear distance R, which we refer to 
as the potential-coupling matrix, and which is written 
in atomic units (a.u.). The atomic unit of energy is the 
Hartree (« 4.36x10“^® J), the atomic unit of length is 
the Bohr radius oq (~ 5.29x10“^^ m), and the atomic unit 
of time is 6 ./i?h (~ 2.42xl0“^^s). The diagonal elements 
of U are the three diabatic potential energies shown in 
Fig. 1 of Lin et al. using the molecular state basis 
given in (IMI) and (IMll . converted to atomic units. The el¬ 
ement is the diabatic radial coupling given in Fig. 3a 
of Ref. [ 83 . The element C /12 is the diabatic rotational 
coupling shown as a dashed line in Fig. 3b of [s^ (and 
incorrectly labeled there as 2 ' ^ 11 + — 1 ' ^H), and C /23 is 
the rotational coupling shown as a solid line in Fig. 3b 
(and incorrectly labeled 1 ' ^ 11 + — 1 ' ^H). 

The potential-coupling matrix is used to define the 
scattering Hamiltonian H, also expressed in atomic units, 
as follows: The diagonal elements are given by 

Ha{R) = Uu{R) + ^{^-^'^. z = 1,2,3 (Cl) 

where Uu^R) are the diabatic electronic potentials dis¬ 
cussed above, /r = 6214.35 a.u. is the Na-He reduced mass, 
b is the impact parameter, and vg is the initial relative 


velocity of the collision (all in atomic units). The addi¬ 
tional centrifugal terms in m are obtained by making 
the classical approximation that the orbital angular mo- 





FIG. 22. (color online) Off-diagaonal elments Ha' of the Na- 
He scattering Hamiltonian with collision parameters vo = 2.0 
a.u. and 6 = 0.5 a.u. 

mentum M is equal to bfivo, with £ ^ 1 and £{£+!) ~ £^. 
Due to the high kinetic energy considered here, the col¬ 
lision dynamics and scattering probabilities are not af¬ 
fected by the centrifugal terms, which only produce an 
energy shift. The diagonal elements of H in SI units are 
plotted in Fig. [21] 

The off-diagonal rotational coupling elements are given 
by 

i/l 2 (R)= X[/I 2 (i?), (C2) 

H23{R) = (^) X U23{R), (C3) 

where the classical approximation to the orbital angular 
momentum is again applied, and the off-diagonal radial 
coupling element is given by 

Hi3iR) = Ui3{R). (C4) 

These are plotted in SI units in Fig. [521 Our definitions of 
U and H follow from Eqs. (4-35) and (4-47b) of Ref. [s^ . 

The i?-dependent scattering Hamiltonian m through 
(IC4I) becomes time-dependent after assuming the semi- 
classical trajectory (l84l) . Scatterers are assumed to have 
initial and final internuclear separations of i? = 50 a.u., 
resulting in the time-dependent Hamiltonian shown in 

Fig. [To! 
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